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We calculate the spectral density of voltage fluctuations in a Single Electron Transistor (SET), 
biased to operate in a transport mode where tunneling events are correlated due to Coulomb interac- 
tion. The whole spectrum from low frequency shot noise to quantum noise at frequencies comparable 
to the SET charging energy (Ec/h) is considered. We discuss the back-action during read-out of a 
charge qubit and conclude that single-shot read-out is possible using the Radio-Frequency SET. 

PACS numbers: 03.67.Lx 42.50.Lc 73.23.Hk 85.25.Na 



The Single Electron Transistor (SET) has recently 
been suggested as a read-out device for solid-state charge 
qubits |l|, ||, 1^, ||] . Presenting new state of the art fig- 
ures for sensitivity in charge measurements, Aassime et 
al. [Q made probable that the Radio-Frequency SET ||] 
(RF-SET) can be used for single shot read-out of the 
Single Cooper-pair Box (SCB) qubit ^. This is pos- 
sible if the measurement time tms needed to resolve the 
two states of the qubit is much shorter than the time 
tmix required to destroy the initial qubit state because of 
back-action due to voltage fluctuations on the SET. 

The mixing rate 1 /tmix is proportional to the noise at 
the frequency w corresponding to the energy splitting be- 
tween the qubit states /S.E. Moreover, in order to have 
well defined charge states to achieve single shot read-out, 
Ai? should be comparable to the qubit charging energy 
Eqb. The ratio Eg^/Ec may be varied between approxi- 
mately 1/10 and 10 in order to optimize the performance 
of the system. 

The voltage fluctuations on the SET island have pre- 
viously been studied in the experimentally accessible low 
frequency limit {hto <^ Ec)- The noise is here fully un- 
derstood both in the classical transport regime of sequen- 
tial tunneling, described by a master equation with fre- 
quency independent tunneling rates |9[ |ll|, as well 
as in the Coulomb blockade cotunneling regime [Q. In 
Ref. g the back-action of the SET was estimated us- 
ing simple interpolation between the low frequency shot 
noise limit and the high frequency limit of Nyquist noise 
from two independent tunnel junctions. 

In this paper we evaluate the finite frequency noise on 
a fully quantum mechanical basis, taking into account 
the frequency dependence of the tunneling rates in the 
presence of fluctuations, separating processes that absorb 
or emit energy. We present a simple expression for the 
voltage noise of a SET biased in the transport mode, 
valid in the whole frequency range from low frequency 
classical shot noise to high frequency quantum noise. By 
comparing tmix with tms for reading out a SCB qubit 
with the RF-SET of Ref. g, we conclude that t^s < 
tmix, i-e. single-shot read-out shoifld indeed be possible. 
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FIG. 1: Schematic figure of the SET capacitively coupled to 
the SCB. 



Consider a small metallic SET island coupled via low 
transparency tunnel barriers to two external leads, and 
coupled capacitively to the SCB, which is controlled by 
a gate voltage (see Fig.|l|). First we calculate the volt- 
age noise on the SET island taking the voltage Vg on 
the SCB to be constant. Then the back-action of these 
fluctuations on the state of the SCB is estimated. This 
approach is appropriate in the considered limit of weak 
SET-qubit coupling {nCc/Cgb, k <Cl; Cl,Cr - Cgb). 

We follow the outline of Ref. @ and model the SET 
by the Hamiltonian 



H = Hl+Hr + Hi + V + Ht = Ho + Ht.. (1) 



where 
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describe noninteracting electrons in the left/right lead 
{Hr,r S {L,R}) and on the island {Hj). The quantum 
numbers n denote transverse channels including spin, and 
k, I denote momenta. The Coulomb interaction on the 
island is described by 

ViN) ^ Ec{N -n, f, (3) 

where N denotes the excess number operator, Ec = 
e^/2C the charging energy (C = Cl + Cr + Cc), 
the fractional number of electrons induced by the exter- 
nal voltages {rix is the fractional part of {ClVl + CrVr + 
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C'cK;)/e) and e the electron charge. The tunnehng term 
is 



N(T) = 11X11 



E E 

r—L,R kin 



(4) 

where the operator e^** changes the excess particle 
number on the island by ±1 and T^" are the tunnel- 
ing matrix elements. $ is the canonical conjugate to 
/V, ([$,iV] = i). In this case of a metallic island, con- 
taining a large number of electrons, the charge degree of 
freedom N = 0, ±1, ... is to a very good approximation 
independent of the electron degrees of freedom l,n. 

The Cooper-pair box is a small superconducting island 
coupled to a superconducting reservoir, via a tunnel junc- 
tion with Josephson energy Ej (see Fig. |l|). The box is 
also characterized by its charging energy Eqb — e^/2Cgb, 
where Cgb is the total capacitance of the box. For 
Ej Eqb two nearby charge states, corresponding to 
zero and one excess Cooper pair on the island, can be 
used as the |0) and |1) states of the qubit. The energy 
splitting AE between these states is controlled by the 
gate voltage Vg. 

For a given charge measurement sensitivity dq and 
measurement time tms , the uncertainty in charge is given 
by Aq — Sq/y^tms- In order to separate the two qubit 
states we need the two intervals ± Aq and 2e ± Aq 
not to overlap, which defines the measuring time tms — 
bq^l{Kef. 

During the measurement, the voltage fluctuations on 
the SET island will induce transitions between the qubit 
states. The rate for excitation/relaxation of the qubit 
is proportional to the SET noise spectral density SyiuS) 
at the negative/positive frequency corresponding to the 
transition AE /%. The information of the initial state is 
destroyed at the combined rate |l^ 



Fi = 
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~ n^"" AE^ 



5.(^) + Sv-(-^) 



, (5) 



here given in the relevant limit AE ^ Ej used below. 

The spectral density of voltage fluctuations on the 
SET island is described by the Fourier transform of the 
voltage-voltage correlation function 



SviLu)^^J dTe-'^^Tr{p,t{to)N{T)NiO)}. (6) 

Here pst{to) is the density matrix of the system in steady- 
state, which is assumed to have been reached at some 
time to before the fluctuation occurs {to < min{0, r}). 
Pst is the tensor product of the equilibrium (Fermi dis- 
tributed) density matrix pl^ for the electron degrees of 
freedom in each reservoir {L, i?, I) and a reduced den- 
sity matrix p^^, describing the charge degrees of freedom, 
which we assume to be diagonal with elements 
denoting the probability of being in charge state N. 
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FIG. 2: a) Diagrammatic expression for S'v(i^) in our approx- 
imation, neglecting processes involving other charge states 
than £ {0,1}. b) Dyson-type equation for Ilff]^i{uj), the 
frequency-dependent transition rate between charge states N 
and N'. 



To evaluate <S'y(ix') we make a perturbation expansion 
of the forward and backward time evolution operators in 
terms of the tunneling term Ht, as shown in Fig. j^a. 
The reservoir degrees of freedom are traced out using 
Wick's theorem. In the diagrammatic language of refer- 
ence (Fig. H) we have forward and backward prop- 
agators of the Keldysh contour (horizontal lines), with 
internal charge transfer vertices e** (small dots) con- 
nected by reservoir lines. The iV(0) and N{t) operators 
form external vertices (big dots). Figure |^a shows the di- 
agrammatic expression for S'y(a;) in our approximation, 
neglecting processes involving other charge states than 
N e {0,1}. Here UNiyf{u}) is the frequency dependent 
rate for transitions between the charge states N and N' . 
One may graphically write down a Dyson type of equa- 
tion for I1nn'{^^) (Fig. ^). In matrix notation this reads 



n(^) 



-n(^) 



muj) 



-E(a>) 



(7) 

where T,n]\ji{uj) represents the sum of irreducible dia- 
grams, i.e. containing no free propagators, for transitions 
between and TV'. 

Since the SET is biased in transport mode, we consider 
only the dominating lowest order diagrams in Y.nj^i(lu), 
corresponding to single tunneling events. In this approx- 
imation we get [E2I 



(8) 



where ^^{10) denotes the rate, in second order pertuba- 
tion theory, to go from charge state to ± 1 while 
the SET is absorbing (uj) or emitting {—lo) a quantum of 
energy \hu}\. The well-known expressions for these rates 
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Ja5n(ac^-A^_i) (10) 



where the dimensionless conductivity 
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PrPl 



Rk 



(11) 



is the ratio of the quantum resistance Rk = h/e^ and 
the resistance of barrier r G {L,R}, and /? = l/ksT. 
n{E) — 1/{1 — exp[— /3£^]} is the Bose-Einstein distribu- 
tion of electron-hole excitations and comes from the con- 
volution of the Fermi distributions for filled initial states 
and empty final states. — V{N)+ iir — V{N — ^ii 
is the energy gained by an electron tunneling from the 
chemical potential of lead r [fir) to the chemical poten- 
tial of the island (/^z), thus changing the charge state 
from to iV + 1 (see Figs. ||a-c). We also assume that 
the relaxation processes on the island are fast on the 
timescales we are looking at and neglect the energy de- 
pendence of both tunneling matrix elements T^" ~ T*"" 
and reservoir densities of states p"{e),p"{e). To avoid 
renormalization effects we choose not too small bias, 
i.e. maxr a5Zn(£'c/|A5|) <C 1. 

We also choose the bias not too high, so that the 
higher charge states are inaccessible at low frequency, 
i.e. 7+(0) = 7o-(0) = giving P^* = for TV ^ {0, 1}. 
At low frequencies we may thus restrict ourselves to the 
two lowest charge states N G {0, 1}. The matrix in- 
verse in Eq. (|^) may then easily be performed analyti- 
cally. For high enough frequencies also the charge states 
N e {—1,2} are accessible. In this regime we may ex- 
pand Eq. (0) to first order in E/w. The results in both 
these regimes combine to the following expression for the 
noise valid, within our approximations, for all frequen- 
cies: 



Sv{uj) 



where PJ* 
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7r(0)/(7o+(0) + 7r(0)) and Pf = 1 
Note the simple structure of this expression: a sum 
over the probabilities of being in the states = and 
iV = 1 times the rates for possible transitions from these 
states, normalized by a denominator containing their fi- 
nite lifetime due to the bias. 

At zero frequency Eq. (p^ ) coincides with the classi- 
cal shot noise result [|[ |10|, p^, which is recovered by 
approximating all rates with their zero frequency value. 
In the high-frequency limit, where the biasing energy is 
negligible, one recovers the high-frequency Nyquist noise 
from two resistances connected in parallell to ground. 
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FIG. 3: a)-c) Schematic illustration of the frequency depen- 
dence of the rate for tunneling from the left lead to the island, 
when the island is in the A*' = state. The shaded area indi- 
cates the amount of electrons energetically allowed to tunnel; 
a) Tiuo = 0, b) Ttw < (SET emits energy) and c) huj > 
(SET absorbs energy), d) Schematic figure showing the bias 
used in the discussion of the noise at different frequencies. 



It is informative to analyze Eq. ( [T^ ) at zero tempera- 
ture where the tunneling rates are linear functions of the 
possible energy gain (since n{E) = Q{E), where 6(x) is 
the unit step function). To be specific we assume that 
PL > Pi > PR and Aq > so that electrons tunnel from 
the left lead to the right and = is the lowest charge 
state. We also assume that the left bias is larger than 
the right one, |A^| > |A^| (see Fig. |a). 

For large negative frequencies hco < - 
tuations try to extract more energy than the bias can 
provide, thus there is no contribution to the noise from 
this frequency range (since we neglect higher order cotun- 
neling effects). Between -|A^| < ^iw < -|A^| the SET 
may emit energy \hLu\ while tunneling through the left 
barrier, while for — |A^| < Tiu < tunneling through 
both barriers contribute. For positive frequencies the 
SET may absorb energy hiu while tunneling through 



|A^| the fluc- 



either barrier. For frequencies tko > |Aq^^| there is 
also a contribution from electrons tunneling backwards, 
against the bias, at the (L/R) junction. In the region 
lAg/^l huj < {\A%^\,\A^^'^\} the higher charge 
states N E {—1,2} are still inaccessible and the noise 
approaches exactly one half of the Nyquist noise. This 
may be explained by the strong correlation of fluctua- 
tions at the two barriers induced by the Coulomb in- 
teraction. At high frequencies huj > {jA^^'^j, |A^^^|} 
also processes at barrier L/R involving the charge states 
N G {—1,2} contribute, giving the real high-frequency 
limit of Nyquist noise from two completely uncorrelated 
tunnel jimctions. Figure ^ illustrates the above discus- 
sion of which processes contribute in which frequency re- 
gion. 

In the low frequency regime \Tiuj\ < |Ag'|, where no 
backward processes are allowed, the noise spectral den- 
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FIG. 4: Schematic picture of the processes contributing to 
the noise in different frequency regimes. Each arrow (passage 
through a barrier) involves emission or absorption of \hu:\. 
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;ity may be written as 

e2 21 /e + 27rtj [P^'a^ + Pfa^] 



C2 



c.2 + 4[7o+(0) + 7r(0)]' 



(13) 



where / = 2e7+(0)7f (0)/(7([(0) +7f (0)) is the DC cur- 
rent through the SET. This is the expression for classical 
shot noise plus a linear w-term in the numerator, indi- 
cating that the probability of processes where the SET 
absorbs/emits energy {±uj) is raised/lowered compared 
to the classical expression (see Figs. ||c,d and the right 
inset of Fig. ^). The total frequency dependent decay 
rate of the states, found in the denominator, is constant 
in this regime since the linear increase / decrease in rates 
for processes absorbing/emitting energy |?ia;| cancel ex- 
actly, at zero temperature. 

The sum of positive and negative frequency SET noise, 
proportional to the back-action on the qubit, is plotted 
in Fig. using the parameters of the RF-SET in Ref. 0| 
with sensitivity 6q = 6.3iie/^/Hz. One may define the 
ignal-to-noise ratio (SNR) of a single shot measurement 
as SNR — y/tmix/tms, which is plotted as a function 
of qubit level splitting AE in the left inset of Fig. ||. 
Using realistic qubit parameters for an SCB made with 
the same aluminum technique as the RF-SET, Eqh/ks — 
IK, AE = 2AEqb, Ej/Eqb — 0.1, assuming read-out 
at AE/Eqb = 2.4 we find that SNR = y^Unix/Uns » 5, 
which indicates that single-shot measurement is possible, 
[f one increases the qubit charging energy, e.g. by using 
a niobium qubit, the SNR increases even further. 

In conclusion, we have obtained a simple expression for 
the finite frequency voltage fluctuations of a SET island 
biased in the transport mode. Using this information 
about back-action we can conclude that single shot read- 
out of a SCB qubit using an RF-SET is possible. 
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5: Summed voltage noise Sv{i^) + Sv{—u) for a symmet- 
biased SET with DC-current 6.7nA and Rt = = 
0.49, T = 20mA:, Ec/kg = 2.57^, Fuh ex- 
;ssion (solid). Classical shot noise (dashed), Nyquist noise 
(dstted). The circles denote, in order from left to right, 
{lAfl = |A«i|, |A«| = |AJ I, |Af I = |A^i|}. The right 
in4et shows the low frequency regime separating positive and 
ne ^ative frequencies (solid) compared with the classical sym- 
metrical expression (dotted). The left inset shows the SNR 
of a single shot read-out as a function of level splitting in the 
qu bit. 
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